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SECTION  1 


INTRODUCTION 

During  the  early-time  evolution  of  high-altitude  nuclear  events, 

1  2 

populations  of  ions  form  with  energies  in  excess  of  10  or  even  100  keV.  ’ 
Subsequently,  a  complex  series  of  collisions  occur  in  which  the  electric 
charge  state  of  the  fast  ions  changes  while  the  ionization  of  the  thermal 
plasma  background  increases.  Simultaneously,  as  their  energy  degrades,  the 
initial  population  of  high  energy  ions  tends  to  become  part  of  the  thermal 
plasma  background.  The  spatial  and  temporal  evolution  of  initial  high  energy 
ionized  structure  must  impact  later  morphology  by  determining,  at  least  in 
part,  the  initial  conditions  for  the  later  fluid  development.  However,  if 
chemical  releases  in  the  ionosphere  are  any  indication,  mechanisms  other  than 
collisions  may  also  influence  the  evolution  of  high  energy  ion  morphology  in 
the  ionosphere.  For  example,  shaped  barium  injection  across  the  magnetic  field 
in  the  ionosphere  has  been  observed  to  cause  the  formation  of  magnetic  field 
aligned  striations  on  an  ion-gyrotime  scale. Called  "prompt  striations" 
in  Reference  3,  the  structures  have  been  attributed  to  a  purely  growing,  finite 
ion-gyroradius,  electrostatic  flute  instability  caused  by  the  ion  loss  cone 
velocity  distribution  function  inherent  to  the  shaped  injection. 

The  primary  purpose  of  this  report  is  to  suggest  that  the  prompt 
striation  phenomena  can  play  a  role  in  the  evolution  of  plasma  jets  crossing 
the  geomagnetic  field  after  a  high  altitude  nuclear  explosion.  In  particular, 
for  a  STARFISH  type  of  nuclear  burst,  prompt  striations  could  develop  at  higher 
altitudes  (i.e.,  well  above  four  hundred  kilometers)  for  superthermal  ions 
with  velocity  less  than  the  maximum  ion  velocity  possible  following  a  nuclear 
burst  explosion.  Another  aim  of  this  report  is  to  point  out  that  a  cyclotron 
resonance  mechanism  can  retard  the  motion  of  both  ions  and  electrons  along 
magnetic  field  lines.  This  cyclotron  resistivity  process  points  out  a  dif¬ 
ference  between  high-altitude  nuclear  and  barium  cloud  environments.  Speci¬ 
fically,  shaped  barium  injection  along  the  ambient  magnetic  field  lines  is 
unaffected  by  the  ion-cyclotron  resistivity  associated  with  Alfven  waves  in 
predominantly  barium  plasmas.  However,  10  to  100  keV  singly  ionized  monatomic 


5 


oxygen  particles  formed  by  high-altitude  nuclear  explosions  and  travelling 
along  the  ambient  magnetic  field  are  susceptible  to  cyclotron  resistivity  for 
reasonable  oxygen  plasma  densities  and  wave  magnetic  field  strength  much 
smaller  than  the  ambient  magnetic  field.  The  implication  of  cyclotron  resis¬ 
tivity  is  that  particle  energy  is  deposited  closer  to  their  source  than  would 
otherwise  be  predicted. 

The  remainder  of  this  report  is  divided  in  the  following  way.  Section 
2  applies  prompt  striation  theory  to  plasma  jetting  across  the  magnetic  field 
under  high-altitude  nuclear  explosion  conditions.  Section  3  describes  the 
cyclotron  resistivity  process  and  a  comparative  application  of  the  process 
to  high-altitude  nuclear  explosion  and  chemical  release  environments.  Section 
4  contains  concluding  remarks  and  points  out  areas  requiring  continued  research. 
Three  self-contained  appendices,  which  support  the  discussion  in  Section  2, 
are  attached  at  the  end  of  this  report.  Specifically,  Appendices  A,  B,  and 
C  discuss  the  effects  of  electron-ion  collisions,  background  plasma  density, 
and  electromagnetic  effects  on  the  formation  of  prompt  striations,  respectively. 


SECTION  2 


PROMPT  STRIATIONS  IN  THE  HIGH-ALTITUDE  NUCLEAR  ENVIRONMENT 

2.1  INTRODUCTION 

Plasma  jetting  across  the  magnetic  field  following  barium  releases 

in  the  ionosphere  has  been  observed  to  cause  the  formation  of  small  scale- 

3  4 

size  field  aligned  structures  on  the  ion-cyclotron  time  scale.  5  These 

structures,  called  prompt  striations  in  Reference  3,  could  impact  the  high- 

altitude  nuclear  problem  in  several  ways.  First,  the  prompt  striations  have 

been  attributed  to  a  kinetic  instability  requiring  finite  ion-gyroradius. 

These  instabilities  are  not  included  in  fluid  models  of  the  high-altitude 

nuclear  environment.  Second,  observation  of  the  prompt  striations  following 

barium  releases  indicates  that  they  are  not  transitory  but  persist  for  periods 

of  time  which  are  much  longer  than  the  time  required  for  their  formation. 

For  example,  in  Avefria  Dos  the  prompt  striations  are  present  for  times  much 

4 

longer  than  one  minute  although  they  are  formed  within  several  seconds.  Third 
the  saturated  state  of  prompt  striations  could  provide  conditions  conducive 
to  subsequent  fluid  growth  once  the  high-velocity  ions,  which  caused  the  ori¬ 
ginal  striation  growth,  become  less  energetic.  Perhaps  this  later  fluid  dyna¬ 
mics  is  the  reason  for  the  prompt  striation  persistence  referred  to  earlier. 

In  this  section  the  prompt  striation  theory  is  applied  to  the  high- 
altitude  nuclear  environment.  Electron-ion  collisions,  cold  background  plasma 
and  finite  ion  dynamic  pressure  all  reduce  the  vigor  of  linear  instability 
growth  and  promote  stability.  These  effects  are  discussed  fully  in  Appendices 
A,  B,  and  C,  and  in  References  5-7. 

2.2  SUMMARY  OF  THE  APPENDICES 

In  this  subsection,  salient  results  found  in  Appendices  A,  B,  and  C 
are  sunmarized.  Particular  attention  is  given  to  those  aspects  of  the  prompt 
striation  theory  which  have  application  to  the  high-altitude  nuclear  explosion 
environment. 

The  linear  instability  postulated  to  be  the  source  of  prompt  striations 
is  a  purely  growing,  finite  ion-gyroradius  flute  mode.  The  free  energy  for 


this  instability  is  supplied  by  high  velocity  ions  crossing  the  ambient  mag¬ 
netic  field.  In  the  analysis  the  ion  velocity  distribution  function  is  assumed 
to  have  the  loss  cone  form: 


^(v)  = 


V±6 


2tr  v/  A(s) 


exp  - 


(Vvo  - 


ir\ 

— ) 


(2.1) 


with 


A ( 6 )  =  (J)*  { (l  +  51 2 3)[1  +  erf (6//2)]  +  (£)*  6  exp(-<52/2)|  , 


(2.2) 


In  Eqs.  (2.1)  and  (2.2),  v,  and  vM  are  the  perpendicular  and  parallel  velocities, 
respectively.  The  parameters  vQ  and  v^  are  constants  which  together  determine 
mean  perpendicular  particle  velocity  and  velocity  spread  about  this  mean  velo¬ 
city.  Because  the  flute  approximation  is  invoked,  the  dispersion  relation  does 
not  depend  on  the  precise  functional  form  for  the  parallel  distribution  func¬ 
tion,  h(v|(),  as  long  as 


oo 

1  ■/  dvu  h<vn> 

— oo 


(2.3) 


For  the  physical  effects  considered  in  the  appendices,  the  dispersion 
relation  and  mode  properties  depend  on  six  dimensionless  parameters: 

1)  6,  which  is  a  measure  of  the  peakedness  of  the  ion  loss  cone 
distribution  function, 

2)  the  linear  growth  rate  in  units  of  the  loss  cone  ion- 
gyrofrequency, 

3)  kvQ/n..,  which  gives  the  wavenumber  perpendicular  to  the  ambient 
field  (k)  in  terms  of  a  typical  ion-gyroradius,  v  , 
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r 
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4)  v  /|n  |,  the  electron-ion  collision  frequency  in  units  of  the 
el ectron-gyrof requency , 

5)  p/p^ ,  the  ratio  of  ambient  to  loss  cone  mass  density,  and 

2  2 

6)  ei  =  (4tt  p .j  vQ  )/B  ,  the  ratio  of  ion  inertial  to  magnetic  field 
pressure  for  loss  cone  ions  with  speed  vQ  and  mass  density  p^. 

Analysis  of  the  dispersion  relation  indicates  that  the  instability 
can  occur  only  in  intervals  of  kvQ/n.  such  that  Jq( k.vQ/!^i )J^(kvQ/ni- )  is  nega¬ 
tive,  where  Jq  and  are  the  ordinary  Bessel  functions  of  the  first  kind  and 
order  zero  and  one,  respectively.  The  interval  of  negative  with  srr  'ist 
kvQ/fii  is  2.5  <,  kvQ/fi^  <,3.8.  The  instability  within  this  interval  of  k 
is  of  greatest  interest  because  it  is  least  susceptible  to  stabilization 
thermal  spread  in  ion  loss  cone  velocity  distribution  function,  electron 
collisions,  and  cold  plasma  background.  The  separate  minimum  electron-u 
collision,  background  plasma  density,  and  criterion  necessary  for  stability 

are  shown  in  Table  2.1  for  6  =  5  and  6  °°.  Values  of  6  =  5  and  6  -*■  °°  are 

specifically  considered  because  the  Buaro  ion  loss  cone  velocity  distribution 
has  been  fitted  to  the  6=5  form  while  the  most  vigorous  instability  occurs 


« 


Table  2.1.  Separate  minimum  electron-ion  collision, 

background  plasma  density,  and  6^  criterion 
necessary  for  stability  are  shown  for 
6  =  5  and  5  -*■  «. 


/  VI".! 

p/p^ 

ei 

6  =  5 

4.2  x  10"3 

.023 

0.2 

6  -*■  °° 

1.5  x  10"2 

.059 

0.56 

?.3  APPLICATION  OF  THE  PROMPT  STRIATION  THEORY  TO  THE  HIGH-ALTITUDE  NUCLEAR 

ENVIRONMENT 

In  real  situations,  a  variety  of  stabilization  mechanisms  can  act 
simultaneously  to  limit  the  growth  of  structure.  The  purpose  of  this  sub¬ 
section  is  to  determine  the  role  of  electron-ion  collision,  cold  plasma  back¬ 
ground,  and  finite  ion  inertia  stabilization  in  relation  to  prompt  striation 
formation  under  high-altitude  nuclear  explosion  conditions.  Microfiche  data 
from  a  Mission  Research  Corporation  simulation  of  STARFISH  is  used  in  the 
discussion.1 

First  the  electron-ion  collision  and  the  ion  inertia  criterion  are 

related  to  one  another  and  then  compared.  Reference  8  shows  that  v  /In  I 

e  1  e‘ 

has  the  following  form: 


with  X  the  Coulomb  logarithm,  me/m..  electron  to  ion  mass  ratio,  c  the  speed 
of  light,  ng  the  density  of  loss  cone  electrons  or  ions,  and  x^  the  electron 
Debye  length.  In  Eq.  (2.4)  ions  have  been  assumed  to  be  singly  charged.  For 

O 

oxygen  ions,  X  a  20,  and  vQ  a  1  x  10  cm/sec,  Eq.  (2.4)  becomes 


.74 


(2.5) 


Now  in  a  plasma  ng  x^  >>  1. 
ng  x ^  >  4.1  x  10 5,  while  if  ng 


6  o 

For  example,  if  n  =  1  x  10  cm  and  Tg  >  1  eV, 
1  x  10®  cm’®  and  Tg  >  1  eV,  ng  x^®  a  4.1 


x  10^.  Moreover,  <  0.56  for  instability.  It  then  follows  from  Eq.  (2.5) 
that  electron-ion  collisional  stabilization  is  unimportant  compared  to  finite 
ion  inertia  stabilization  under  early-time  high  altitude  nuclear  explosion 
conditions.  Even  if  vQ  is  1  x  10®  cm/sec  rather  than  1  x  10®  cm/sec,  finite 
ion  inertia  stabilization  would  generally  take  precedence. 


O 

Consider  singly  ionized  oxygen  ions  with  vQ  =  1  x  10  cm/sec.  From 
Table  2.1  it  follows  that  the  ion  inertia  constraint  permits  instability  for 
oxygen  loss  cone  ion  densities  satisfying: 
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10 


6  =  5 


(2.6a) 


6  -*■  «  .  (2.6b) 

In  Eqs.  (2.6a)  and  (2.6b),  the  magnetic  field  strength  is  in  units  of  gauss. 
Also,  Eqs.  (2.6a)  and  (2.6b)  and  Table  2.1  imply  that  instability  requires 
cold  background  oxygen  ion  densities  satisfying: 

B2  cm'3  ,  6=5  (2.7a) 

B2  cm"3  ,  6  +  -  .  (2.7b) 

Figures  2.1  -  2.8  are  taken  from  microfiche  provided  by  Reference  1 
and  represent  relevant  output  from  the  simulation  of  STARFISH.  The  com¬ 
bination  of  these  figures  and  Eqs.  (2.6)  and  (2.7)  suggests  that  high  velocity, 

O 

vQ»  1  x  10  cm/ sec,  ions  are  not  likely  to  be  the  source  of  prompt  striations 

in  high-altitude  nuclear  events  like  STARFISH  for  three  reasons.  First,  at 

low  altitude  both  the  plasma  density  and  magnetic  field  strength  are  largest. 

However,  the  maximum  magnetic  field  strength  in  the  blast  wave  is  roughly  1 

gauss  while  the  electron  number  density  at  these  locations  is  greater  than 
5  -3 

1  x  10  cm  .  Second,  at  higher  altitudes  the  magnetic  field  strength  is 

3  -3 

approximately  0.3  gauss  while  the  electron  density  exceeds  3  x  10  cm  where 

o 

ion  velocity  is  1  x  10  cm/sec.  If  cold  plasma  number  density  is  greater  than 
2-3  2 

1.3  x  10  cm  for  6  =  5  or  9  x  10  for  6  -*■  ®,  instability  is  impossible.  Such 
cold  plasma  number  densities  seem  likely.  Third,  at  approximately  two  seconds 

Q 

mean  ion  velocity  is  less  than  1  x  10  cm/sec  (see  Figure  2.8). 

To  circumvent  the  finite  ion  inertia  stabilization  which  is  so  effec¬ 
tive  when  the  magnetic  field  is  weak  and  the  ion  velocity  is  large,  it  is 
necessary  to  consider  whether  prompt  striation  formation  might  be  easier  if 
the  ion  velocity  is  reduced.  For  example,  if  vQ  =  1  x  10^  cm/sec  for  singly 
ionized  oxygen,  then  finite  ion  inertia  and  cold  plasma  stabilization  imply 
that 
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Figure  2.3.  Electron  number  density  (cm-3)  at  150.2  milliseconds 
for  the  STARFISH  VER238  simulation 
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gure  2.5.  Magnetic  field  line  geometry  at  2000.7  milliseconds 
for  the  STARFISH  VER238  simulation 
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Figure  2.6.  Magnetic  field  strength  (gauss)  at  2000.7  milliseconds 
for  the  STARFISH  VER238  simulation 
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Figure  2.7.  Electron  number  density  (cm"3)  at  2000.7  milliseconds 
for  the  STARFISH  VER238  simulation 
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Figure  2.8.  Mean  ion  velocity  (cm/sec)  at  2000.7  milliseconds 
for  the  STARFISH  VER238  simulation 
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Clearly,  Eqs.  (2.8)  and  (2.9)  are  more  easily  satisfied  than  Eqs.  (2.6)  and  (2.7) 

g 

appropriate  to  vq  *  1  x  10  cm/sec  loss  cone  ions.  Furthermore,  Figures  2.3 
and  2.7  show  that  plasma  density  generally  decreases  with  increasing  altitudes. 
Hence,  prompt  striations  are  not  likely  to  develop  from  the  highest  energy 
ions  at  low  altitude,  below  400  kilometers,  but  rather  from  ions  with  more 
modest  energies  at  higher  altitudes,  well  above  400  kilometers. 


SECTION  3 


ION  AND  ELECTRON  CYCLOTRON  RESISTIVITY 

3.1  INTRODUCTION 

In  space  and  laboratory  plasmas,  electric  fields  and  wave  particle 
processes  can  cause  the  formation  of  a  superthermal  electron  or  ion  population 
whose  momentum  is  aligned  primarily  along  the  ambient  magnetic  field  lines. 
Despite  the  inability  of  classical  resistivity  to  counter  the  acceleration  of 
these  fast  particles  under  many  circumstances,  plasma  mechanisms  exist  which 
tend  to  impede  acceleration  along  magnetic  field  lines  to  infinite  energies. 
Landau  damping  of  bump-in-the-tail  instabilities  is  an  example  of  such  an  ano- 

g 

malous  resistive  plasma  process.  Also,  in  strongly  magnetized  plasmas  with 
the  electron-cyclotron  frequency  smaller  than  the  electron  plasma  frequency, 
the  anomalous  electron-cyclotron  resonance  can  markedly  inhibit  the  motion 
of  fast  electrons  along  the  ambient  magnetic  field. 

The  cyclotron  resonance  process  is  shown  to  be  an  effective  anomalous 
resistivity  mechanism  for  ions  as  well  as  electrons  in  both  strongly  and  weakly 
magnetized  plasmas  if  two  conditions  are  satisfied.  First,  it  is  necessary 
that  the  wave  frequency  be  smaller  than  the  relativistically  correct  gyrofre- 
quency.  In  fact,  strong  anomalous  resistivity  can  occur  for  waves  with  fre¬ 
quency  approaching  zero  although  energy  transfer  between  waves  and  particles 
occurs  on  a  much  slower  time  scale.  Second,  the  wave  electric  and  magnetic 
fields  perpendicular  to  the  ambient  magnetic  field  must  have  an  appreciable 
wavenumber  parallel  to  the  ambient  field  which  satisfies  the  cyclotron  res¬ 
onance  condition.  An  implication  of  the  theory  is  that  cyclotron  resistivity 
requires  shorter  parallel  wavelengths  for  strongly  magnetized  plasmas  than  for 
weakly  magnetized  plasmas. 

Examples  of  cyclotron  resistivity  are  presented  which  emphasize  ion¬ 
ospheric  applications;  however,  the  theory  also  has  wider  application  to  space 
and  laboratory  plasmas.  Shaped  barium  injection  along  the  ambient  magnetic 
field  is  demonstrated  to  be  unaffected  by  the  ion-cyclotron  resistivity  process 

with  Alfven  waves  in  predominantly  barium  plasmas  for  reasonable  barium  den- 

3  13-14 

sities  in  the  ionosphere.  ’  However,  10  to  100  keV  singly  ionized  monatomic 
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oxygen  particles  travelling  along  the  ambient  magnetic  field  are  shown  to  be 
susceptible  to  cyclotron  resistivity  for  reasonable  oxygen  plasma  densities 
and  wave  magnetic  field  strength  much  smaller  than  the  ambient  magnetic  field. 
Such  high  energy  ions  can  result  from  high-altitude  nuclear  events  *  as  well 
as  geophysical  processes.15,16 


3.2  THEORETICAL  MODEL  FOR  CYCLOTRON  RESISTIVITY 

To  demonstrate  the  cyclotron  resistivity  process,  we  use  an  infinite, 
homogeneous,  relativistical ly  correct  quasil inear  diffusion  model.  In  parti¬ 
cular,17 
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In  Eq.  (3.1)  the  particle  distribution,  f,  is  a  function  of  the  momentum  com¬ 
ponents  perpendicular  to  and  parallel  to  the  umbient  magnetic  field  (i.e., 
p^  and  p((,  respectively).  In  Eqs.  (3.1)  and  (3.2)  subscripts  x,  y,  and  II 
denote  vector  components  in  the  x,  y,  and  magnetic  field  directions,  respec¬ 
tively,  and  the  symbol  e  represents  a  unit  vector.  The  ambient  magnetic  field 
strength,  particle  charge,  particle  mass,  velocity,  and  speed  of  light  are 
given  by  B,  q,  m,  v,  and  c,  respectively.  For  electrons  or  negative  ions, 
the  charge  is  negative.  The  wavenumbers,  k,  and  the  frequency,  w,  are  related 
through  the  dispersion  relation.  The  summation  in  the  index  k  is  over  modes 
and  Jn  is  the  Bessel  function  of  the  first  kind  and  order  n. 

Because  this  analysis  emphasizes  cyclotron  resonance  effects  for  par¬ 
ticles  with  momentum  which  is  primarily  magnetic  field  aligned,  the  particle 
distribution  function  is  assumed  to  have  the  following  form: 

f(pi>,,,i)  =  8(pi>  (3-3) 

for  real  space  density  nQ.  In  Eq.  (3.3),  6(p  )  and  h(p(|)  are  the  Dirac  delta 
function  in  p^  and  a  general  function  of  p|(,  respectively. 

With  the  choice  of  particle  distribution  function  given  by  Eq.  (3.3), 
it  follows  that  only  the  n  =  -I,  0  and  +1  cyclotron  resonances  can  be  effec¬ 
tive  in  altering  particle  trajectories  because  they  do  not  require  finite 
gyroradius.  The  subsequent  discussion  will  emphasize  the  n  =  -1  and  +1 
resonances  which,  unlike  the  n  =  0  Landau  resonance,  can  be  very  weakly 
dependent  on  wave  frequency. 

To  demonstrate  that  the  n  =  -1  and  +1  resonances  can  be  effective 
anomalous  resistivity  mechanisms,  it  is  sufficient  to  take  appropriate  moments 
of  Eq.  (3.1).  For  an  arbitrary  function  of  p^  and  p|(,  g(p^,p|(),  the  rate  at 
which  the  g ( p ^ » P , ( )  moment  changes  for  mode  k  and  cyclotron  harmonic  n  is: 
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In  Eqs.  (3.5b)  and  (3.5d), 
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which  is  derivable  from  Eq.  (3.1)  and  Faraday's  law.  In  Eq.  (3.6),  wave 
magnetic  fields  are  denoted  by  B„  and  B  . 


Equation  (3.5d)  shows  that  wave  fields  must  vary  with  time  if  energy 
is  to  be  transferred  between  fields  and  charged  particles.  For  nonzero  fre¬ 
quencies,  Eqs.  (3.5)  and  (3.6)  also  indicate  that  the  n  =  +1  resonance  causes 
energy  to  be  removed  from  negatively  charged  particles  and  to  be  added  to  posi¬ 
tively  charged  particles.  Similarly,  the  n  =  -1  resonance  adds  and  subtracts 
energy  from  negatively  and  positively  charged  particles,  respectively. 

Irrespective  of  how  or  whether  energy  is  exchanged  between  electro¬ 
magnetic  or  electrostatic  fields  and  particles,  Eqs.  (3.5)  and  (3.6)  demon¬ 
strate  that  the  parallel  momentum  is  always  reduced  while  the  perpendicular 
momentum  becomes  greater  for  both  the  n  =  -1  and  +1  resonances  as  long  as 
to  <  |£2/y|.  The  key  requirement  for  pitch  angle  scattering  is  the  presence 
of  particles  satisfying  the  resonance  condition, 

k  p  /m  =  yw  -  nfi  (3.7) 

or  for  yw  <<  |nfl | , 

k(|p|)/m  as  -  nfl  .  (3.8) 

Equation  (3.8)  shows  that  as  particle  parallel  momentum  increases,  the  parallel 
wavelength  required  for  cyclotron  resistivity  becomes  longer. 


3.3  COMPARATIVE  APPLICATION  BETWEEN  NUCLEAR  AND  BARIUM  CLOUD  ENVIRONMENTS 

The  purpose  of  this  subsection  is  to  illustrate  a  difference  between 
high-altitude  nuclear  and  barium  cloud  environments  in  the  ionosphere.  Speci¬ 
fically,  in  the  example  presented  here,  it  is  demonstrated  that  10  to  100  keV 
oxygen  ions  injected  parallel  to  the  magnetic  field  lines  by  nuclear  explosion 
can  be  retarded  by  ion-cyclotron  resistivity  with  Alfven  waves  for  reasonable 
densities,  e.g.,  1  x  10  cm'  ;  however,  barium  ions  injected  parallel  to  the 
magnetic  field  by  shaped  injection  are  not  retarded  by  the  ion-cyclotron 
resistivity  with  Alfven  waves  in  a  predominantly  barium  plasma  unless  the 
density  is  relatively  high  (i.e.,  much  greater  than  1.4  x  10  cm  ).  Such 
high  barium  densities  may  be  difficult  to  achieve  for  the  shaped  injection. 


For  the  cyclotron  resistivity  process  to  be  effective  for  waves  with 
frequencies  well  below  the  ion-cyclotron  frequency,  Eq.  (3.8)  must  be  satis¬ 
fied.  Figure  3.1  shows,  for  B  =  0.3  gauss,  the  necessary  wavelengths  for 
resonance  as  a  function  of  barium  and  monatomic  oxygen  energy.  The  figure 
has  three  significant  features.  First,  parallel  wavelength  increases  with 
ion  energy.  Second,  for  equal  energies  barium  requires  longer  parallel  wave¬ 
lengths  to  satisfy  Eq.  (3.8)  than  oxygen.  Third,  the  parallel  wavelengths 
necessary  for  cyclotron  resistivity  are  in  excess  of  1  kilometer  for  ion 
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energies  greater  than  100  eV.  For  1  x  10  eV  ions,  the  required  parallel 

2 

wavelengths  are  approximately  1  x  10  kilometers  for  barium  and  40  kilometers 
for  oxygen. 

As  a  specific  quantitative  example  we  consider  the  case  of  Alfven 
waves  propagating  along  the  ambient  magnetic  field  and  interacting  with  the 
superthermal  ions  travelling  along  the  ambient  magnetic  field.  The  origin 
of  the  Alfven  waves  is  not  important  for  the  present  example  but  could  be 
associated  with  a  nuclear  explosion,  barium  release,  or  other  abrupt  change 

13  1  c 

in  the  ionosphere.  ’  Besides  having  a  frequency  below  the  ion-cyclotron 
frequency,  the  significant  feature  of  the  waves  is  that  for  propagation  parallel 
to  the  ambient  magnetic  field,  the  wave  has  no  parallel  electric  field  associated 
with  it.  Hence,  only  the  wave  magnetic  field  contributes  to  cyclotron  resis¬ 
tivity  in  Eq.  (3.5). 

Consider  a  barium  ion  with  a  kinetic  energy  of  160  eV  and  speed  of 
1.5  x  10^  cm/sec.  For  B  =  0.3  gauss,  Eq.  (3.8)  and  Figure  3.1  indicate  that 
a  parallel  wavelength  of  4.5  kilometers  is  required  for  cyclotron  resistivity 
to  alter  the  ion  pitch  angle.  For  the  wave  frequency  to  be  much  smaller  than 
the  barium-cyclotron  frequency  it  is  necessary  that 

k„vA  «  ni  (3.9) 

with  vA  the  Alfven  speed  and  n.  the  ion-cyclotron  frequency.  With  B  =  0.3 

gauss,  Eq.  (3.9)  implies  that  the  resonance  condition  cannot  be  satisfied 

7  -3 

unless  the  barium  ion  density  is  much  greater  than  1.4  x  10  cm  .  Hence, 
unless  such  high  barium  densities  can  be  made  to  exist  for  shaped  barium 
injection  along  the  magnetic  field  in  the  ionosphere,  cyclotron  resonance 
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The  solution  of  Eq.  (3.8),  parallel  wavelength 
in  kilometers  versus  ion  energy  in  electron  volts, 

for  barium  ( - )  and  monatomic  oxygen  ( - ). 

The  magnetic  field  strength  is  0.3  gauss. 


between  Alfven  waves  and  superthermal  barium  ions  in  a  predominantly  barium 
plasma  is  highly  improbable. 

4  5 

Monatomic  oxygen  ions  with  kinetic  energy  of  1  x  10  and  1  x  10  eV 

can  undergo  cyclotron  resistivity  with  B  =  0.3  eV  only  if  2tt/ k )(  =  12  and  38 

kilometers,  respectively.  The  monatomic  oxygen  densities  required  to  satisfy 

Eq.  (3.9)  for  the  aforementioned  parallel  wavelengths  are  in  excess  of  2.3 
“>  4-3 

x  10“  and  2.3  x  10  cm  ,  respectively.  Hence,  in  contrast  to  the  barium 
example  presented  before,  oxygen  cyclotron  resistivity  due  Alfven  waves  is 
possible  in  the  ionosphere  if  the  wave  magnetic  field  strength  is  large 
enough. 

When  it  is  active,  cyclotron  resistivity  is  significant  only  if  the 
momentum  drag  caused  by  it  is  larger  than  classical  collisional  drag.  For 
monatomic  oxygen  energies  greater  than  3  keV  and  an  electron  temperature  of 
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0.1  eV,  Coulomb  drag  on  electrons  is  much  larger  than  Coulomb  drag  on  ions. 

Consequently,  for  cyclotron  resistivity  due  to  mode  k  and  cyclotron  harmonic 

18 

n  to  dominate  over  classical  resistivity  it  is  necessary  that 

^  'Vnh.n  =•  5'4  x  10'55  Vo  f  %  P»'2 

•  00 

h(p(|)  ln(Aie)  ,  (3.10) 

_3 

where  ng  is  the  electron  density  (cm  )  and  the  Coulomb  integral,  ln(A.  ), 
is  specified  later.  The  right-hand  side  of  Eq.  (3.10)  is  valid  so  long  as 
the  distribution  function,  h(p(|),  corresponds  to  most  oxygen  particles  having 
energies  greater  than  3  keV. 

To  permit  quantitative  evaluation  of  the  left-hand  side  of  Eq.  (3.10), 
we  assume  the  following  form  for  h ( p (( ) : 

h(p„)  =  (2ir  m.2v.2)  i  exp[-(p)(  -  P0)2/2mi2vi2]  (3.11) 

with  pQ  and  v^  constants  and  m^  the  ion  mass. 
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Conventional  plasma  instability  theory  might  suggest  that  the  speci¬ 
fied  form  of  h(Pj|)  could  be  susceptible  to  waves  which  grow  via  the  inverse 
Landau  damping  mechanism,  especially  if  the  parallel  phase  velocities  of  the 
waves  is  less  than  p0/m.j .  We  circumvent  this  possibility  by  noting  that  for 
the  chosen  distribution  function  inverse  Landau  damping  will  not  promote  in¬ 
stability  for  parallel  phase  velocities  in  excess  of  pQ/m. .  In  any  case, 
cyclotron  resistivity  can  act  whether  or  not  another  source  of  anomalous 
resistivity  is  active. 

The  right-hand  side  of  Eq.  (3.11)  can  be  integrated  approximately  by 
the  saddle  point  methods  if  the  distribution  function  is  strongly  peaked  with 
po  >;>  mivi  and  the  weak  dependence  of  the  Coulomb  logarithm  on  momentum  is 
neglected.19  Equation  (3.10)  becomes: 

(Viih.n  >  5'4  *  10'55  Vo  V2 

exp(-fQ)  ln(A.e)  (3.12) 


with 


fQ  =  0.5  52  (yx  -  l)2  +  2  ln(yi)  ,  fQ' '  =  62  -  2/y  *  , 

=  °-5  tl  +  (1  -  S/S2)*]  ,  5  =  p0/m.v.  .  (3.13) 

Writing  the  argument  of  the  Coulomb  logarithm  as:20 


permits  the  complete  specification  of  the  right-hand  side  of  Eq.  (3.13).  In 
Eq.  (3.14),  the  parameters  e,  k,  Tg,  m..,  and  mg  denote  the  magnitude  of  the 
electron  charge,  Boltzmann's  constant  electron  temperature,  ion  mass,  and 
electron  mass,  respectively. 
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Figure  3.2  is  a  plot  of  wave  magnetic  field  strength  in  gauss  versus 

2 

monatomic  oxygen  energy  (0.5  pQ  /m. )  in  eV  necessary  to  satisfy  Eq.  (3.12) 
for  n  =  1  x  106  cm’^,  B  =  0.3  gauss,  T  =  0.1  eV,  6=5,  and  m.n./k  p  =  1.5. 

C  v  C  II  II  v 

The  figure  shows  that  wave  magnetic  field  strengths  much  smaller  than  the  am¬ 
bient  magnetic  field  are  sufficient  to  permit  cyclotron  resistivity  to  dominate 

over  classical  resistivity  especially  for  high  energy  oxygen  ions.  In  particu- 
2 

lar,  for  p  /2m.  values  of  10  keV  and  100  keV,  wave  magnetic  field  strengths 

01  -3-4 

for  comparable  cyclotron  and  classical  resistivity  are  5  x  10  and  9  x  10 

gauss,  respectively.  Larger  values  of  the  wave  fields  result  in  the  dominance 

of  cyclotron  resistivity. 


Plot  of  fluctuation  magnetic  field  strength 
i By  ±  iBx |  ;c  B  (gauss)  versus  oxygen  kinetic 
energy  for  B  =  0.3  gauss,  ne  =  1  x  105  cnr3, 
and  Te  =  0.1  eV.  Shaded  region  above  the  curve 
indicates  the  parameter  space  for  which  cyclotron 
resistivity  dominates  over  classical  resistivity. 


SECTION  4 


CONCLUDING  REMARKS 


Prompt  striations  have  been  observed  in  some  barium  release  experi- 
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ments  as  plasma  tries  to  jet  across  the  ambient  magnetic  field.  ’  They  have 

been  postulated  to  be  the  result  of  a  collective  process  driven  by  free  energy 

available  from  an  ion  loss  cone  distribution  function  associated  with  the 
3  5-7 

jetting.  ’  It  has  been  shown  in  the  appendices  that  electron-ion  colli¬ 
sions,  background  plasma  density,  and  finite  ion  inertial  pressure  all  tend 
to  reduce  the  linear  growth  rate  of  the  instability.  However,  for  application 
of  the  theory  to  the  high-altitude  nuclear  environment  electron-ion  collisional 
stabilization  is  unimportant.  Finite  ion  inertial  pressure  seems  to  virtually 
eliminate  the  possibility  that  the  most  energetic  ions  formed  after  a  high- 
altitude  explosion  can  be  the  source  of  prompt  striations  in  the  weak  magnetic 
field  environment  of  the  ionosphere.  However,  in  a  burst  like  STARFISH,  prompt 

striations  could  develop  from  a  class  of  ions  with  lower  velocities  perpendi- 
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cular  to  the  ambient  magnetic  field,  e.g.,  between  1  x  10  and  1  x  10  cm/sec. 

The  development  of  prompt  striations  would  benefit  from  the  lower  background 
plasma  density  environment  available  at  higher  altitudes,  especially  well  above 
400  kilometers  for  a  STARFISH  type  burst. 

Ion  and  electrons  moving  along  magnetic  field  lines  have  been  shown 
to  be  susceptible  to  a  pitch  angle  scattering  process  associated  with  a  cyclo¬ 
tron  resonance  phenomena.  An  implication  of  this  c'clotron  resistivity  is 
that  high  energy  ions  and  electrons  travel  less  rapidly  along  magnetic  field 
lines  than  would  otherwise  be  predicted  and  so  must  deposit  their  energy  closer 
to  their  source.  It  has  been  shown  that  the  particular  process  of  ions  inter¬ 
acting  with  Alfven  waves  can  more  easily  occur  in  the  high-altitude  nuclear 
rather  than  in  the  barium  cloud  environment.  Hence,  ion-cyclotron  resistivity 
with  Alfven  waves  represents  a  way  that  barium  releases  can  differ  from  nuclear 
explosions  in  the  ionosphere. 

Questions  still  remain  concerning  the  applicability  of  prompt  striation 
and  cyclotron  resistivity  processes  to  the  high-altitude  nuclear  environment. 
Regarding  prompt  striations,  the  kinetic  nonlinear  saturated  state  is  not  obvious 
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in  the  ionosphere.  An  interesting  possibility  is  that  the  linear  instability 
postulated  to  be  the  source  of  prompt  striations  readily  couples  to  longer 
wavelengths  in  guiding  center  fashion.  In  Reference  21  it  is  suggested  that 
the  energy  spectrum  of  the  system  might  go  as  k  for  long  wavelengths  and  k~ 
for  short  wavelengths.  This  guiding  center  saturation  mechanism  should  be 
compared  to  other  possible  saturation  processes  for  the  instability.  It  is 
also  important  to  know  whether  the  correct  kinetic  saturated  state  could  form 
the  basis  for  later  fluid  evolution  as  ion  energies  decrease  and  fluid 
assumptions  become  more  reasonable.  Regarding  cyclotron  resistivity,  the 
ultimate  momentum  distribution  function  which  results  from  interaction  between 
waves  and  particles  has  not  been  determined.  This  distribution  function  is 
necessary  for  determining  the  spatial  and  temporal  energy  deposition  by  super- 
thermal  particles  by  other  processes,  say  charge  exchange.  The  possibility 
that  field  aligned  structuring  could  develop  as  pitch  angle  scattering  occurs 
should  also  be  studied. 
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APPENDIX  A 


ELECTRON-ION  COLLISION  EFFECTS  ON  PROMPT  STRIATION  FORMATION 

ABSTRACT 

Electron-ion  collisions  have  been  included  in  the  analysis  of  electro¬ 
static,  purely  growing,  finite  ion-gyroradius,  flute  instabilities  generated 
by  a  loss  cone  distribution  of  energetic  ions  in  the  ionosphere.  It  is  shown 
that  the  collisions  not  only  reduce  the  growth  rate  below  the  value  predicted 
by  collisionless  theory,  but  also  result  in  a  second  lower  growth  rate  mode 
not  found  in  collisionless  theory.  Although  the  electron-ion  collision  fre¬ 
quency  is  larger  than  the  ion-gyrofrequency  for  the  parameters  of  the  Buaro 
barium  release  experiment,  the  growth  rate  of  the  linearly  most  unstable  mode 
is  not  appreciably  modified  by  the  collisions  from  the  value  predicted  from 
collisionless  theory.  Electron-ion  collisions  have  a  greater  impact  on  the 
instability  as  the  electron  density  is  increased.  Complete  stabilization 
through  electron-ion  collisions  requires  densities  in  excess  of  8.7  x  107 


A-l.  INTRODUCTION 


The  formation  of  magnetic  field  aligned  structures,  striations,  is 

an  ionospheric  physics  problem  of  considerable  interest  because  of  their 

A1  A2 

adverse  impact  on  communication.  At  low  altitudes  a  generally  agreed 

upon  process  for  striation  growth  is  the  EXB  gradient  drift  instability  (e.g., 

see  Reference  A3  for  a  detailed  analysis).  However,  at  higher  altitudes 

striations  are  not  readily  attributed  to  the  gradient  drift  mechanism  but 

A4 

are  nevertheless  observed  in  shaped  barium  releases  such  as  Buaro.  The 
authors  of  Reference  A4  call  the  structures  observed  in  Buaro  prompt  stria¬ 
tions  because  of  their  rapid  appearance  and  attribute  these  striations  to  an 
electrostatic  finite  ion-gyroradius,  purely  growing,  flute  instability  generated 
by  an  ion  loss  cone  distribution  function  appropriate  for  the  plasma  following 
the  barium  release.  Corroboration  for  the  theory  of  Reference  A4  is  the  experi¬ 
mental  observation  in  Buaro  of  various  ion-cyclotron  harmonic  modes. 

The  purpose  of  this  paper  is  to  include  electron-ion  collisions  into 
the  analysis  of  the  finite  ion-gyroradius  flute  instability.  Motivation  for 
considering  electron-ion  collision  effects  is  that  even  for  the  electron  den¬ 
sities  in  the  Buaro  experiment^4  (i.e.,  4  x  10^  cm~^),  the  electron-ion  colli¬ 
sion  frequency  is  much  larger  than  the  ion-gyrofrequency  which  is  the  charac¬ 
teristic  growth  rate  for  the  finite  ion-gyroradius,  purely  growing,  modes. 

It  is  shown  that  electron-ion  collisions  do  not  markedly  alter  the  instability 
dispersion  relation  fc  Buaro  parameters  but  could  have  a  significant  effect 
for  higher  density  experiments. 

This  appendix  is  divided  into  four  sections.  In  Section  A-2  the 
theoretical  model  and  some  qualitative  implications  of  the  model  are  described. 
In  Section  A-3  quantitative  solutions  of  the  dispersion  relation  are  presented 
and  discussed.  Section  A-4  is  a  discuss, on  of  the  significant  features  of 
the  analysis  and  results. 


A-2.  THEORETICAL  MODEL 


The  theoretical  model  assumes  an  infinite  homogeneous  plasma,  with 
a  loss  cone  distribution  of  ions  superposed  on  a  Maxwellian  distribution  of 
electrons.  Reference  A4  indicates  that  this  is  the  appropriate  distribution 
function  after  the  barium  release.  For  the  purely  growing  modes  being  con¬ 
sidered  the  inverse  wavenumber  ts  smaller  than  the  ion-cyclotron  radius  but 
larger  than  the  electron-cyclotron  radius.  Consequently,  the  ion  response 
to  the  modes  requires  a  fully  kinetic  treatment  while  the  electron  response 
can  be  recovered  from  the  fluid  equations.  Electrostatic  flute  modes  with 
the  plane  wave  variation,  exp[i(k^y  -  wt}],  are  also  assumed  with  k^,  y,  w, 
and  t  being  the  wavenumber  component  perpendicular  to  the  ambient  magnetic 
field,  y-directed  coordinate,  wave  frequency  in  radians  per  second,  and  time, 
respectively.  With  the  assumed  plasma  and  wave  model,  the  ion  response  is 
identical  to  that  found  in  Reference  A4.  Specifically: 


bi 


E 

P=1 


2u 


(2p-l)  I 


2  n  2o  2 

u)  "  P  C- 


p.p-i 


2  ,  \2„  '2 
U>  ”  (p-l  ) 


(A.  1) 


where 


!P,P-1 
A  ( 6 )  » 


? 

477  n^e  /mi  ,  ^  =  eB^uuc  ,  bi  =  vo/fi.  , 

=  jjjj  f  d*  x  Jp(k1bix)  Jp.1(kibix)  exp[-0. 5<52(x-l )2] 

0 

(7)^(1  +  62)[1  +  erfCd/v7?)]  +  (|)*6  exp(-0.562)}  (A. 2) 


In  Eq.  (A. 2),  the  parameters  n..,  e,  m. ,  BQ,  and  c  denote  the  ion  density, 
magnitude  of  the  elementary  electron  charge,  ion  mass,  magnetic  field  strength, 
and  speed  of  light,  respectively.  The  ordinary  Bessel  function  of  the  first 


kind  and  order  p  is  represented  by  J  .  For  the  singly  ionized  quasineutral 
plasma  assumed  in  the  present  calculation, 
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The  parameters  vQ,  v..,  and  5  =  vQ/v^  are  associated  with  the  form  for  the  loss 
cone  ion  distribution  function, 
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where  v^  and  v|(  are  the  perpendicular  and  parallel  velocity  components, 
respectively. 


The  electron  response  to  the  electrostatic  flute  mode  is: 
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In  the  expression  for  the  electron  collision  frequency,  which  is  an  approxima¬ 
tion  valid  near  thermal  equilibrium,  the  density  and  temperature  are  in  units 
of  cm”3  and  eV,  respectively/3  Equation  (A. 6)  includes  both  the  electron 
polarization  and  Pederson  current  responses  to  the  wave  electric  field  and 

can  be  derived  either  from  kinetic  equations  (e.g.,  Eq.  11  of  Reference  A6 
2  2  2 

in  the  limit  of  k.  v  /fi  -+  0  where  v„  is  the  electron  thermal  speed)  or  from 
X  s  e  e 
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fluid  equations.  An  assumption  made  in  the  present  analysis  and  justified 
in  Subappendix  AA  is  that  the  ion  response  to  waves  is  unaffected  by  electron- 
ion  collisions  while  the  electron  dynamics  in  the  wave  are  not  modified  by 
the  ion  response. 

The  dispersion  relation  can  now  be  written  as: 


For  purely  growing  modes,  iy  =  w  with  y  real.  With  the  definitions, 


y  *  y/ni  ,  a  =  vQ/u. 


(A. 8) 


and  the  general  ionospheric  relation. 
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the  dispersion  relation  can  be  rewritten  in  the  following  way: 
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For  a  •>  0,  this  expression  is  equivalent  to  the  homogeneous  plasma  dispersion 
relation  found  in  Reference  A4. 

For  m../me  >>  1,  a  good  approximation  to  Eq.  (A. 10)  is: 
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where  |n  |  is  the  magnitude  of  the  electron-cyclotron  frequency.  The  signi¬ 
ficant  feature  of  Eq.  (A. 11)  is  that  for  specified  values  of  vg,  kb.,  and  6, 
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all  ion  species  have  the  same  eigenvalue,  y.  However,  the  absolute  growth 
rate  is  larger  for  less  massive,  higher  gyrofrequency  ion  species.  For 
example,  the  growth  rate  for  monatomic  oxygen  is  greater  than  for  barium  by 
the  ratio  of  the  barium  to  oxygen  particle  masses. 
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A-3.  RESULTS 


In  the  absence  of  electron-ion  collisions  and  a  very  peaked  distri¬ 
bution  function  with  vQ  »  v^,  there  are  a  very  large,  but  finite  denumerable, 
number  of  zones  in  k^b..  space  which  permit  positive  growth  rate  solutions  of 
the  dispersion  relation.  The  mathematical  reason  for  these  different  zones, 
which  correspond  to  negative  values  of  J ^ ( k^b ^ ) ,  is  the  sinusoidal 

nature  of  Bessel  functions  and  the  slow  decrease  in  maximum  amplitude  of  the 
Bessel  function  oscillations  as  k^b^  increases.  The  maximum  possible  nor¬ 
malized  growth  rate  [(y/q. )m,  1  versus  k.b.  is  shown  for  the  first  three 
unstable  zones  in  Figure  A.l  (zone  I,  2.5  <  k,b.  <  3.8;  zone  II,  5.6  <  k,b.  < 

7.0;  zone  III,  8.7  <  k.b^  <_  10.1).  A  characteristic  feature  of  each  of  the 
solutions  is  the  sharp  increase  in  growth  rate  near  the  marginal  values  of 
k^.  which  permit  growth.  The  change  in  growth  rate  is  not  as  noticeable  for 
values  of  kJtK  away  from  the  marginal  values.  The  figure  also  shows  that  the 
maximum  growth  rate  of  each  zone  increases  with  zone  number,  at  least  for  the 
first  three  zones  illustrated. 

In  the  subsequent  discussion  detailed  quantitative  results  will  be 
confined  to  the  first  two  unstable  zones.  Additional  qualitative  comments 
will  be  made  concerning  the  higher  number  unstable  zones. 

Equation  (A. 10)  has  been  evaluated  numerically  to  determine  the  electron- 
ion  collision  frequency  consistent  with  linear  wave  growth  rate  and  perpendicu¬ 
lar  wavenumber.  Figures  A. 2a  to  A.2f  are  plots  of  normalized  collision  frequency 
versus  normalized  linear  wave  growth  rate  for  various  values  of  kjb..  and  5. 

A  general  characteristic  feature  of  each  curve  in  the  figures  is  the 
two  nonzero  values  for  the  growth  rate  which  can  occur  when  the  electron-ion 
collision  frequency  is  not  equal  to  zero.  In  particular,  for  zero  collision 
frequency  the  only  solution  of  Eq.  (A. 10)  which  yields  a  positive  growth  rate 
is  the  one  determined  by  the  dispersion  relation  of  Reference  A4.  As  the  elec¬ 
tron-ion  collision  frequency  increases,  the  maximum  growth  rate  of  the  instability 
decreases  below  the  value  calculated  neglecting  collisions.  In  addition,  a 
second  purely  growing  solution  of  lower  growth  rate  exists.  This  solution 
is  not  found  in  a  theoretical  model  which  neglects  collisions.  As  the  electron- 
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ible  normalized  growth  rate,  (Y/fi-j)max>  versus 
°°  and  the  first  three  unstable  zones 
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ion  collision  frequency  is  increased  further,  the  growth  rates  of  the  two  solu¬ 
tions  to  Eq.  (A. 10)  become  more  comparable.  Eventually,  for  the  maximum  colli¬ 
sion  frequency  which  permits  solution  to  Eq.  (A. 10)  for  a  specific  kb.,  the 
two  solutions  have  the  same  growth  rate.  This  growth  rate  is  lower  than  the 
maximum  growth  rate  calculated  in  the  absence  of  collisions,  but  nevertheless 
can  be  a  substantial  fraction  of  the  ion-cyclotron  frequency  (i.e.,  y/r^.  ,>0.1). 

Comparing  the  plots  in  Figures  A. 2a  to  A.2f,  it  is  clear  that  as  the 
parameter  5  gets  larger  the  maximum  collision  frequency  which  permits  solution 
to  Eq.  (A. 10)  increases  and  the  instability  gets  more  vigorous.  The  physical 
mechanism  behind  this  increased  instability  strength  is  the  greater  availa¬ 
bility  of  velocity  space  free  energy  as  the  velocity  space  distribution  function 
becomes  more  singular.  Mathematically,  the  magnitude  and  sign  of  the  Bessel 
functions  for  a  specific  value  of  kb..  is  not  washed  out  by  integration  over 
the  distribution  function. 

Figures  A. 2a  to  A.2f  demonstrate  that  the  linear  instability  growth 
and  the  variation  of  growth  rate  with  the  electron-ion  collision  frequency 
both  depend  in  a  sensitive  way  on  the  magnitude  of  .  For  example,  in 
Figure  A. 2a  the  maximum  growth  rates  for  k^b..  *  3.1  and  3.6  are  approximately 
0.55  and  0.3  respectively.  The  k^b^  *  3.1  mode  is  stabilized  by  Vg/S^. 

>,  980,  while  the  kj.b.  =  3.6  mode  is  stabilized  by  much  lower  values  of  the 
collision  frequency,  Vg/fi^  >,87. 

Figures  A. 3a  and  A. 3b  show  the  maximum  value  of  Vg/fl..  [i.e.,  (ve/n1-)max3 
which  permits  a  solution  of  Eq.  (A. 10)  for  the  first  two  unstable  zones  and 
specified  value  for  5.  Values  of  Vg/n,.  which  exceed  the  magnitude  indicated 
by  the  curves  cannot  satisfy  Eq.  (A. 10),  while  those  values  below  the  curve 
can  satisfy  the  equation.  Unstable  zone  I  solutions  have  been  found  for 
<5  >  4  but  not  for  5  _<  3;  however,  for  zone  II,  6  _>  6  provides  a  solution 
while  5  <_  5  does  not.  Since  Reference  A4  primarily  considered  6=5,  that 
paper  does  not  include  the  second  unstable  zone. 

Comparison  of  Figures  A. 3a  and  A. 3b  indicates  that,  although  Figure  A.l 

shows  that  the  zone  II  instability  can,  in  principle,  have  a  greater  growth 

rate  than  the  zone  I  instability,  electron-ion  collisions  stabilize  the  zone 

II  instability  more  easily  than  the  zone  I  instability.  For  6  ”,  (v./fl. )_.w 

s  i  max 

is  approximately  equal  to  3750  in  zone  I  and  1100  in  zone  II.  Qualitative 
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Figure  A. 3a.  The  maximum  value  of  Vg/n.. ,  (ve/ni)max>  which  permits 
solution  of  Eq.  (A. 9)  versus  6  for  zone  I.  (ve/fi-j)max 
for  <s  -+  »  is  indicated  on  the  right-hand  side  of  the 
plot.  Values  of  Vg/ft^  which  exceed  the  magnitude  indi 
cated  by  the  curves  cannot  satisfy  Eq.  (A. 10),  while 
those  values  below  the  curve  can  satisfy  the  equation. 
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Figure  A. 3b.  The  maximum  value  of  v./a,,  (vfl/Ji.)ma  .  which  permits 

si  ci  max 

solution  of  Eq.  (A. 9)  versus  <5  for  zone  II.  (v  /a.) 
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for  6  -*-»  is  indicated  on  the  right-hand  side  of  the 
plot.  Values  of  ve/ft..  which  exceed  the  magnitude  indi¬ 
cated  by  the  curves  cannot  satisfy  Eq.  (A. 10),  while  those 
values  below  the  curve  can  satisfy  the  equation. 
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examination  of  the  third  and  higher  unstable  zones  shows  the  general  trend 
that  stability  through  electron-ion  collisions  is  more  easily  achieved  as 
k  b.,  increases.  The  physical  reasoning  behind  this  stabilization  behavior 
is  that  waves  which  require  higher  harmonics  of  the  wavelength  to  fit  within 
the  circumference  of  an  ion-cyclotron  orbit  are  more  easily  disrupted  by 
phenomena  which  prevent  ion-cyclotron  orbits  from  being  perfect  circles. 

stabilization  for  large  k^b^  is  easier  because  charac- 
function  oscillations  about  zero  decrease  in  magnitude 
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A-4.  DISCUSSION 


By  including  electron-ion  collisions  in  the  analysis  of  electrostatic, 
finite  ion-gyroradius,  purely  growing,  flute  modes,  it  has  been  demonstrated 
that,  when  unstable  solutions  to  the  dispersion  relation  exist,  there  are 
generally  two  solutions  which  satisfy  the  dispersion  relation.  Both  solu¬ 
tions  have  growth  rates  which  lie  below  the  vaiues  predicted  by  collisionless 
theory.  Above  the  value  of  the  collision  frequency  for  which  both  solutions 
merge,  no  purely  unstable  mode  can  exist  for  the  specified  values  of  k^ 
and  6. 

A4 

The  paper  by  Simons,  et  al.,  proposes  that  the  finite  gyroradius, 
purely  growing  mode  is  the  source  of  prompt  striations  in  the  Buaro  barium 
release  and  so  it  is  useful  to  determine  the  role  of  electron-ion  collisions 

with  the  assumed  and  measured  parameters  of  the  Buaro  plasma  (i.e.,  n.  as 

5-3  * 

4  x  10  cm"  ,  B0  «  0.32  gauss,  and  6  *  5).  For  an  electron  temperature  of 

Tga  0.1,  the  ratio  of  vg/n.  is  equal  to  21.7.  Examining  Figure  A. 2a,  it  is 
evident  that  the  electron-ion  collisions  in  the  Buaro  plasma  do  not  signifi¬ 
cantly  alter  the  growth  rate  of  the  most  unstable  mode  predicted  by  collision- 
less  theory  (i.e.,  k^.  a  3.1  and  y/n^  w  0.55). 

For  plasma  densities  which  are  higher  than  appropriate  to  the  Buaro 

experiment,  the  electron-ion  collision  frequency  has  a  corresponding  greater 

stabilizing  impact  on  the  finite  gyroradius  purely  growing  flute  modes.  For 

7  -3 

example,  if  n^  as  lx  10  cm  and  Tg  »  0.1  eV,  then  vg/Q^  as  480  and  the  maxi¬ 
mum  growth  rate  for  k^  as  3.1  in  Figure  A. 2a  is  reduced  to  0.48  from  the 
collisionless  value  of  0.55  .  Figures  A. 3a  and  A. 3b  show  that  with  ve/«i  a: 

480,  no  solution  is  possible  for  6  <,  4  in  zone  I  and  6  <,  11  in  zone  II.  For 

7-3 

even  higher  densities,  n^  ,>  8.7  x  10  cm  ,  Figures  A. 3a  and  A. 3b  show  that 
self  electron-ion  collisions  are  sufficient  to  completely  stabilize  the  finite 
gyroradius,  purely  growing,  modes.  Such  relatively  high  electron  densities 
are  required  to  stabilize  the  instability  because  of  the  weak  electron  polari¬ 
zation  and  Pederson  response  to  the  modes. 
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SUBAPPENDIX  AA 


In  this  subappendix  it  is  shown  that  the  lowest  order  ion  response 
to  the  waves  is  unaffected  by  electron-ion  collisions,  while  the  electron 
dynamics  in  the  wave  are  not  modified  by  the  ion  response. 

The  demonstration  that  the  lowest  order  ion  response  is  not  affected 
by  electron-ion  collision  consists  of  two  parts.  First,  the  characteristic 
frequency  for  ions  to  change  their  momentum  by  a  substantial  amount  due  to 
collisions  with  electrons  is: 

v-j  «  (me/m^ )  ve  •  (AA.l) 

For  typical  plasma  densities  and  temperatures  in  the  ionosphere  (i.e.,  na  < 
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10  cm  and  Tg  «  0.1  eV),  the  quantitative  value  of  is  very  much  smaller 

than  the  characteristic  linear  growth  rate  of  the  instability  (i.e.,  y  ,<  s^.). 

Second,  the  electron  response  to  the  waves  [see  Eq.  (A. 5)]  is  in  magnitude 

very  much  smaller  than  the  individual  terms  in  the  ion  response  [see  Eq.  ( A . 1) ] . 

Hence,  the  electron  currents  in  the  waves  cannot  significantly  alter  the  ion 

dynamics  through  electron-ion  collisions. 

To  show  that  the  ion  coupling  through  electron-ion  collisions  does 
not  significantly  modify  the  electron  response  to  the  waves  under  considera¬ 
tion,  it  is  convenient  to  examine  the  electron  momentum  equation. 


-ii#  n  v  . 
e  e  el 


=  -eEj  ne 


-  m  n  v  (v  , 

e  e  e'  el 


-«u> 


(AA.2) 


where  subscript  "1"  denotes  wave  quantities.  To  lowest  order  in  inverse  mag¬ 
netic  field  strength. 
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To  next  order, 


(i)  * 

iijjn  n  v  .  =  en  ~  1-- - -  +  m  n  v  (v  ^  -  v-, )  .  (AA.4) 
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From  Eq.  (AA.4)  it  is  evident  that  the  ion  velocity  in  the  waves  does  not 
significantly  affect  the  electron  velocity  if: 


!vil/vel 


«  1  . 


(AA.5) 


Now  the  dielectric  elements  which  reflect  the  ExB  velocity,  in  the  limit  of 
low  frequency  and  infinitesimal  wave  number,  are  <x  ^  and  k x,  where  x  and 
y  are  the  two  coordinates  orthogonal  to  the  magnetic  field.  In  particular, 
for  electrons: 


iKx.y(6)l  '  'V.x<e)l  =  l>2/Vl  • 


(AA.6) 


while  for  the  ions: 
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For  purely  growing  modes  iy  =  oi  and  it  follows  that: 
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with  the  equality  being  valid  in  the  limit  of  y  -*■  0.  In  deriving  Eq.  (AA.8), 
the  standard  Bessel  function  addition  theorem. 


1  *  J„2(z)  +  2  £  Jp2(z) 


(AA.9) 


P=1 
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has  been  used. 


In  the  limit  of  k  -►  0,  Eq.  (AA.8)  becomes: 
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Examining  Eqs.  (AA.6)  and  (AA.9),  it  is  evident  that: 
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with  the  equality  holding  if  y.k^  -►  0. 

For  the  values  of  frequency  and  wavenumber  considered  in  this  paper 
(i.e.,  growth  rates  a  sizeable  fraction  of  the  ion-cyclotron  frequency  and 
k  v  /«i  >  2.5),  the  inequality  holds  and  Eq.  (AA.5)  is  satisfied  to  lowest 
order. 
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STABILIZATION  OF  PROMPT  STRIATIONS  BY  BACKGROUND  PLASMA 


ABSTRACT 


The  ambient  background  plasma  is  shown  to  have  an  important  stabi¬ 
lizing  influence  on  the  generation  of  electrostatic,  purely  growing,  finite 
ion-gyroradius,  flute  instabilities  driven  by  ion  loss  cone  distribution 
functions  associated  with  shaped  barium  injection  across  the  ambient  magnetic 
field  in  the  ionosphere.  This  stabilization  implies  that  such  short  scale 
size  instabilities  and  the  striations  they  produce  will  only  be  present  if 
the  background  plasma  density  is  low  enough.  Consistent  with  the  experi¬ 
mental  observations  of  the  striations,  it  is  shown  that  for  the  parameters 
of  the  Buaro  barium  release  experiment  the  background  plasma  density  is 
insufficient  to  completely  stabilize  the  finite  ion-gyroradius  modes;  however, 
background  ion  density  does  reduce  the  growth  rate  from  the  level  predicted 
with  a  theory  which  does  not  include  the  background  plasma  dynamics. 


B-l.  INTRODUCTION 


Plasmas  containing  a  population  of  ions  with  a  loss  cone  velocity  dis¬ 
tribution  function  are  susceptible  to  instabilities  which  drive  ions  into  the 
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loss  cone  at  a  much  faster  rate  than  classical  collisional  processes.  ’ 

In  fact,  these  instabilities  are  thought  to  be  the  limiting  factor  for  modern 
B3-B5 

mirror  machines.  “  However,  loss  cone  instabilities  are  not  isolated  to 
the  laboratory. 

Recently  it  has  been  suggested  that  electrostatic,  purely  growing, 
finite  ion-gyroradius,  flute  instabilities  caused  by  an  ion  loss  cone  distri¬ 
bution  function  are  responsible  for  the  prompt  striations  observed  in  the  Buaro 

R6 

shaped  barium  release.  The  purpose  of  this  paper  is  to  show  that  these  finite 
ion-gyroradius  instabilities  are  not  a  general  feature  of  ionospheric  plasma 
expansion  across  a  magnetic  field.  In  fact,  they  can  be  stabilized  for  reali¬ 
zable  values  of  the  background  plasma  density  even  though  complete  stabilization 
does  not  occur  for  the  parameters  of  the  Buaro  experiment.  Quantitative  results 

presented  here  are  not  only  consistent  with  those  of  Swift  for  the  special 

B7 

case  of  monoenergetic  loss  cone  ions  but  also  include  the  effects  of  back¬ 
ground  plasma  for  more  distributed  loss  cone  ion  distribution  functions. 
Moreover,  they  are  in  agreement  with  research  in  the  mirror  fusion  program 
which  indicates  enhanced  stability  of  loss  cone  modes  when  less  energetic 
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plasma  is  present.  ’  When  stabilization  does  occur,  instabilities  char¬ 
acterized  by  slower  growth  times  and  larger  scale  lengths  such  as  the  EXB 
gradient  drift  instability  [e.g.,  see  Reference  B8  for  a  detailed  analysis] 
must  determine  striation  structure  rather  than  the  characteristically  more 
vigorous  and  smaller  scale  size  finite  ion-gyroradius  modes. 

This  appendix  is  divided  into  four  sections.  Section  B-2  describes 
the  theoretical  model  for  the  calculation  and  qualitative  implications  of  the 
dispersion  relation.  Section  B-3  contains  quantitative  results  of  the  solu¬ 
tion  to  the  dispersion  relation,  while  Section  B-4  summarizes  the  theory  and 
results.  Suggested  experiments  for  testing  the  theory  are  presented. 


B-2.  THEORETICAL  MODEL 


An  infinite  homogeneous  plasma  model  consisting  of  two  components  is 
assumed.  The  first  component  is  the  thermal  plasma  background  which  corres¬ 
ponds  to  the  ambient  ionosphere.  The  second  component  is  the  high  velocity 

ion  component  (with  associated  electrons)  which  arises  from  the  shaped  barium 
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release  or  possibly  some  other  natural  or  artificial  phenomena,  ’  and  is 

assumed  to  have  a  loss  cone  velocity  distribution  shown  to  be  the  appropriate 

B6 

distribution  function  after  a  shaped  barium  release. 

The  thermal  background  plasma  consists  of  electrons  with  density  ngo 
and  several  singly  ionized  ion  species  (designated  by  subscript  j)  of  density, 
n.  .  Quasineutral ity  requires, 

J  * 

%=  E  %  •  o-1' 

j 


The  background  magnetic  field  has  strength  BQ.  The  modes  considered  in  this 
appendix  are  plane  waves  corresponding  to  the  electrostatic  limit  of  extra¬ 
ordinary  mode  propagation.  They  are  assumed  to  be  purely  growing  with  growth 
rate  y,  and  to  be  flute-like  with  the  only  spatial  variation  perpendicular 
to  the  magnetic  field  (k.).  The  electric  potential  associated  with  the  wave 
satisfies  Laplace's  equation,  k^  [1  +  4*  (<^  +  Kg)]^  *  0,  where  and 
are  the  ambient  and  loss  cone  plasma  susceptibilities,  respectively. 

With  the  chosen  plasma  model  the  contribution  of  the  ambient  plasma 
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to  the  electrostatic  dielectric  is: 
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where 


“eo2  -  4*  neo  e2/”e  •  “j02  *  4*  "jo  e2/mj  >  bj  *  Yaj  • 


Vj  “  *  ne  =  ‘eBo/n,ec  »  aj  =  eB0/mjc  ‘  (B.3) 


In  Eqs.  (B.2)-(B.3),  the  symbol  I  represents  the  modified  Bessel  function 
of  the  first  kind  and  order  p.  Also,  e,  m  ,  m-,  T.,  and  c  represent  the  mag- 
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nitude  of  the  electron  charge,  the  electron  mass,  mass  of  a  single  particle 
of  background  ion  species  j,  temperature  of  background  ion  species  j,  and  the 
speed  of  light,  respectively.  Two  assumptions  made  in  writing  Eq.  (B.2)  are 
that  the  electron-gyroradius  is  much  smaller  than  the  inverse  perpendicular 
wavenumber  and  that  the  wave  frequency  is  much  smaller  than  the  electron- 
gyrofrequency.  Also,  collision  effects  are  neglected. 

The  high  velocity  ion  plasma  contribution  to  the  dielectric  is  given 
by;B6,B10 

k2  3  (4ir)-1 

where 

we  “  4ir  n.  e2/me  ,  u,.2  =  m.u>e2/me  , 
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In  Eq.  (8.5)  the  parameters  and  m..  denote  the  density  of  loss  cone  ion  species 

i  and  the  mass  of  a  single  particle  of  loss  cone  ion  species  i,  respectively. 

The  ordinary  Bessel  function  of  the  first  kind  and  order  p  is  represented  by 

J  .  The  parameters  v  .  v.,  and  6  =  v„/v,.  are  associated  with  the  assumed  form 
p  0  1  0  1 

for  the  loss  cone  ion  distribution  function, 
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Here  and  v)(  are  the  perpendicular  and  parallel  velocities  to  the  magnetic 
field,  respectively.  The  precise  functional  form  of  g(v())  is  unimportant  for 
the  flute  modes  being  considered. 

The  wave  electric  field  potential  is  nonzero  only  if  the  dispersion 
relation  is  satisfied: 


0  =  1  +  4tt  (k^  +  Kg)  •  (B.7) 

The  new  feature  of  this  dispersion  relation,  which  has  not  been  considered 
previously  for  the  barium  cloud  environment,  is  the  contribution  of  the 
ambient  plasma. 

The  significance  of  the  background  plasma  upon  the  linear  dynamics 
of  the  purely  growing  finite  ion-gyroradius  flute  modes  arises  because  such 
modes  tend  to  be  quasineutral  for  ionospheric  parameters  and  the  ion  dynamics 
in  response  to  the  wave  fields  are  much  more  profound  than  the  weaker  electron 
polarization  currents.  Quantitatively, 
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Hence,  a  very  good  approximation  to  the  dispersion  relation  is: 
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Examination  of  Eq.  (B.9)  demonstrates  that  the  background  ions  always 
promote  stability  because  these  ions  introduce  a  positive  definite  contribution 
to  the  dispersion  relation  while  the  fast  ion  portion  of  the  dispersion  relation 
can  be  positive  or  negative  but  is  always  finite.  Indeed,  the  numerical  eval¬ 
uation  of  Eq.  (B.9)  also  indicates  that  for  a  specified  value  of  k^  and  loss 
cone  ion  density,  growth  rate  is  reduced  as  background  ion  density  is  increased. 
Hence,  lack  of  solution  at  marginal  stability, 
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is  sufficient  to  assure  no  growing  solutions  are  possible  for  a  specified 
k^bi  and  background  plasma  density.  In  deriving  Eq.  (B.IO),  the  identity 
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has  been  used  while 
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A  final  simplification  to  the  dispersion  relation  can  be  made  by  noting 
that  the  weighted  ambient  background  mass  density,  F,  reduces  to  the  actual 
ambient  mass  density  when  b  .  -*•  0  (i.e.,  F  bT^Q  \  n.  m.).  This  aPProx imation 
can  be  very  reasonable  for  high  velocity  ions  (particularly  heavy  ions  like 
barium)  within  the  ambient  ionosphere  because 

bJ./bi  «  1  ;  (B .  13) 

however,  care  must  always  be  taken  in  applying  the  bj  0  approximation  since 
positive  growth  rate  solutions  to  Eq.  (B.9)  are  characterized  by: 


2.5  <  kJLbi  .  (B.14) 

A  general  statement,  deduced  from  Eqs.  (B.10)  and  (B.12),  is  that  maxi¬ 
mum  stabilization  by  the  ambient  plasma  takes  place  when  the  background  is 
sufficiently  cold  so  that  the  kb.  «  1  approximation  is  valid  and  F  has  its 
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largest  possible  magnitude  for  given  loss  cone  and  background  ion  density. 
Physically,  this  most  effective  stabilization  arises  because  the  background 
plasma  current  in  the  potential  fluctuations  is  largest  when  k^  -*•  0  and  is 
optimally  able  to  short-circuit  the  tendency  for  wave  potential  growth  asso¬ 
ciated  with  the  loss  cone  ions.  However,  when  the  thermal  spread  in  the  back¬ 
ground  ion  distribution  function  is  large,  background  ion  current  is  reduced 
and  higher  background  density  is  required  for  stability.  Ultimately,  stability 
for  k^bj  large  is  caused  by  the  reduction  of  the  net  destabilizing  positive 
slope  for  the  total  (i.e.,  loss  cone  plus  background)  ion  distribution  function. 
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B-3.  RESULTS 


Results  from  the  quantitative  evaluation  of  the  dispersion  relation, 

Eq.  (8.9),  are  presented  in  Table  B.l  and  Figures  B.l  and  B.2.  As  noted  in 
Reference  89,  there  are  under  ideal  conditions  with  6  «,  a  very  large  but 

finite  number  of  possible  unstable  zones  in  kjb.  which  correspond  to  negative 
J^k^b^ )Jo(k^bi).  In  the  first  column  of  Table  B.l,  the  zones  with  k^.  <  26 
are  listed.  For  each  zone  of  instability  in  k^b^,  the  value  of  k.b^  which 
is  most  difficult  to  stabilize  with  cold  background  ions  is  given  in  the 
second  column  and  the  corresponding  minimum  value  of  p/p.  required  for  sta¬ 
bility  is  shown  in  the  third  column. 

The  key  features  of  Table  B.l  are  that  instability  is  not  possible 
for  p/p.  >  0.059  and  that  smaller  values  of  p/p.  are  required  for  stability 
as  the  value  of  kb.  in  column  2  increases.  These  results  are  consistent  with 
those  of  Reference  B7  for  a  monoenergetic  ion  population.  Mathematically, 
this  reduction  comes  about  because  the  maximum  amplitude  of  Bessel  function 
oscillations  decreases  with  argument.  Physically,  the  polarization  drift  of 
the  background  ions  in  the  waves  most  easily  stabilize  fluctuations  requiring 
the  highest  degree  of  coherence  between  ions  and  waves  (i.e.,  waves  with  large 
k^bj).  A  corollary  remark  is  that  instability  with  sizeable  k^  requires  6 
large  in  order  that  the  ion  motion  within  the  wave  has  the  necessary  coherence. 

Curves  for  the  marginally  stable  values  of  p/p^  versus  kjb.  are  plotted 
in  Figures  B.la  -  B.lc  for  6=5,  10,  and  «,  respectively.  These  cases  are 
not  quantitatively  analyzed  in  Reference  B7.  In  Figure  B.la  only  the  2.5  < 
k  b^  <  3.8  instability  is  shown  because  for  6  =  5  it  is  the  only  zone  in 
k  b|  which  permits  instability.  For  6  =  10  and  °°,  other  unstable  zones  in 
k  b^  are  possible,  but  for  clarity  only  the  first  two  are  illustrated  in 
Figures  B.lb  and  6.1c. 

The  curves  in  Figures  B.la  -  B.lc  characteristically  have  a  maximum 
magnitude  and  decrease  on  both  sides  of  this  point.  The  rate  of  decrease  is 
most  rapid  for  k^  near  the  values  for  which  JQ(k  )  ^(kjb.)  =  0.  Only 
p/p.  below  the  curves  permit  instability.  In  particular,  as  the  magnitude 
of  p/p.  increases,  only  those  k.b.  which  are  not  too  close  to  the  J  (k.b.) 
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Figure  B.la. 


6  =  5 


7/p^  versus  for  5=5.  The  illustrated 
curve  gives  the  marginal  stability  condition. 

For  values  of  p/p^  below  the  curve,  instability 
is  possible,  but  above  the  curve  modes  are  stable. 


Figure  B.lb.  p/p.,  versus  for  6  =  10.  The  illustrated 


curves  give  the  marginal  stability  condition. 


For  values  of  p/p.  below  the  curves,  instability 
is  possible,  but  above  the  curves  modes  are  stable. 


Figure  8.1c.  p/p^  versus  k  for  6  ->  ».  The  illustrated 
curves  give  the  marginal  stability  condition. 

For  values  of  p/p^  below  the  curves,  instability 
•  is  possible,  but  above  the  curves  modes  are  stable. 
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J 1 ( kb i )  =  0  points  can  be  unstable.  For  example  in  Figure  B-lb,  p/p.  -  0.03 
implies  that  instability  can  occur  only  if  2.6  <,  k^b ^  <  3.4.  For  values  of 
J/o.  above  the  curves  the  criterion  given  by  Eq.  (B.10)  is  satisfied  and  there 
is  stability.  The  three  figures  demonstrate  that  instability  is  possible  for 
larger  values  of  7/Pi  if  <5  increases. 

In  Reference  B6  a  barium  cloud  with  6=5  was  considered.  Figure  B.2 
is  a  plot  of  p/p.  versus  y/^  for  the  corresponding  most  unstable  purely  grow¬ 
ing  mode  in  the  absence  of  background  ion  density  (i.e.,  y/n.  »  0.55  and 
k  b.j  «  3.1).  Given  a  value  for  p/p..,  the  corresponding  solution  to  the  dis¬ 
persion  relation  is  found  at  the  intersection  with  the  plotted  curve.  It  is 
evident  that  with  a  finite  background  ion  density  the  growth  rate  is  reduced 
below  what  it  would  be  if  p/p.-  ■*  0.  Assuming  a  background  electron  density 
of  ng0  =  3  x  10  cm"  ,  a  barium  ion  density  of  n.  =  4  x  10  cm,  and  singly 

ionized  monatomic  oxygen,  then  y  «  0.39  n. .  Although  this  value  for  the  growth 
rate  is  still  a  sizeable  fraction  of  the  ion-gyrofrequency,  the  background 
density  has  resulted  in  a  30%  reduction  in  growth  rate.  Complete  stabilization 

by  background  ion  density  is  not  achieved  in  Buaro.  Figure  B.2  indicates  that 
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a  background  monatomic  oxygen  density  in  excess  of  6.9  x  10  cm  would  be 

required  to  totally  stabilize  the  kjb^  =  3.1  and  6  =  5  instability.  From 

Figure  B-la  total  stabilization  for  6=5  and  all  kb.  is  possible  for  mona- 
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tomic  oxygen  densities  larger  than  8  x  10  cm  .  This  density  is  less  than 
the  barium  density  and  readily  occurs  in  the  ionosphere. 

The  theoretical  model  and  the  application  to  Buaro  have  neglected  mag¬ 
netic  field  strength  and  plasma  spatial  variations.  Reference  B6  shows  that 
Buaro  density  gradients  perpendicular  to  the  magnetic  field  only  weakly  modify 
the  infinite  and  homogeneous  analysis  of  the  purely  growing  loss  cone  insta¬ 
bility,  but  no  comments  are  made  about  effects  associated  with  inhomogeneities 
parallel  to  the  background  magnetic  field.  Now,  a  local  approximation  for 
mode  characteristics  is  reasonable  if  the  plasma  permits  modes  with  very  long 

wavelengths  parallel  to  the  magnetic  field.  In  particular,  the  dispersion 
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relation,  Eq.  (B.2),  is  asymptotically  valid  if: 
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where  x()  is  the  parallel  wavelength.  For  Buaro  parameters,  Eq.  (B .  15 )  implies 
xil  >:>  °*8  km*  A1so«  Figure  3  of  Reference  B6,  it  is  clear  that  the  longi¬ 

tudinal  dimensions  of  the  barium  cloud  are  well  in  excess  of  10  kilometers. 
Hence,  a  calculation  using  local  values  of  density  and  magnitude  field  strength 
is  justified  to  lowest  order. 
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B-4.  SUMMARY  AND  DISCUSSION 


Purely  growing,  finite  ion-gyroradius,  flute  instabilities  generated 
because  of  an  ion  loss  cone  distribution  function  can  be  stabilized  by  the 
ions  of  a  sufficiently  high  density  ambient  ionosphere.  Since  the  electron 
polarization  currents  in  the  waves  are  much  smaller  than  the  ion  currents, 
electrons  do  not  contribute  in  a  significant  way  to  the  wave  dynamics. 

The  purely  growing  finite  ion-gyroradius  flute  instabilities  have 
characteristic  scale  lengths  and  growth  times  which  are  much  smaller  than 
those  typical  of  EXB  gradient  drift  instabilities.  Hence,  the  finite  ion- 
gyroradius  instabilities  provide  a  lower  limit  on  striation  scale  length  in 
the  ionosphere  than  can  be  predicted  from  EXB  gradient  drift  theory  and  can 
be  the  dominant  mechanism  producing  striations  when  operative  (i.e.,  with  the 
background  plasma  density  low  enough). 

B6 

For  the  parameters  of  the  Buaro  shaped  barium  release,  the  local 
approximation  for  the  instability  has  been  justified  and  it  has  been  demon¬ 
strated  that  the  instability  exists  albeit  with  a  linear  growth  rate  reduced 
by  the  background  plasma.  The  theory  with  background  plasma  dynamics  included 

does  indicate  that  if  the  ambient  ion  density  (assumed  to  be  primarily  mona- 
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tomic  oxygen)  had  been  larger  than  8  x  10  cm  rather  than  the  measured  value 
4  -3 

of  3  x  10  cm  ,  the  finite  ion-gyroradius  flute  modes  would  not  have  been 
unstable  and  small  scale  length  prompt  striations  produced  by  these  modes 
would  not  have  been  observed. 

The  theory  and  results  presented  in  this  paper  suggest  two  experiments 
to  test  the  theory.  First,  shaped  barium  releases  similar  to  Buaro  but  with 
higher  background  ion  density  (e.g.,  greater  than  8  x  10  cm  for  monatomic 
oxygen)  should  not  result  in  striations  with  scale  lengths  characteristic  of 
the  finite  ion-gyroradius  flute  instabilities.  Second,  if  a  shaped  barium 
release  is  followed  in  approximately  five  seconds  by  one  which  is  not  shaped, 
the  theory  predicts  that  the  striations  should  first  appear,  if  the  background 
density  is  low  enough,  but  should  eventually  disappear  in  conjunction  with 
ionization  of  the  second  barium  release.  Specifically,  the  Buaro  striations 
would  have  been  completely  stabilized  by  modest  values  of  thermal  barium  ion 
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density  (i.e.,  greater  than  9  x  10  cm  ).  Since  for  both  experiments  the 
finite  ion-gyroradius  instability  is  or  becomes  inoperative,  ionosphere 
striation  structure  formed  must  "freeze11  at  a  larger  scale  length  typical 
of  more  fluid-like  modes  such  as  the  EXB  gradient  drift  instability. 
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APPENDIX  C 


ELECTROMAGNETIC  EFFECTS  IN  THE  ANALYSIS  OF  PROMPT  STRIATIONS 


ABSTRACT 


In  a  number  of  barium  injection  experiments  in  the  ionosphere,  field 
aligned  structures  have  been  observed  to  develop  on  the  ion-cyclotron  time 
scale.  These  structures,  called  prompt  striations,  have  been  attributed  to 
plasma  jetting  across  the  ambient  magnetic  field  and  have  previously  been 
analyzed  in  the  electrostatic  limit.  The  electrostatic  analysis  is  reasonable 
for  experiments  in  which  the  directed  ion  kinetic  energy  per  unit  volume  is 
much  smaller  than  the  magnetic  field  pressure;  however,  for  higher  altitude 
experiments  in  which  the  ambient  magnetic  field  strength  is  weaker,  justi¬ 
fication  for  the  electrostatic  treatment  fails  and  the  full  electromagnetic 
dispersion  relation  must  be  examined.  It  is  demonstrated  that  electromagnetic 
effects  inhibit  the  growth  of  the  linear  instability  considered  to  be  the  source 
of  prompt  striations  in  the  earlier  experiments.  Consistent  with  the  electro¬ 
static  analysis,  cold  background  plasma  and  thermal  spread  in  the  ion  loss 
cone  distribution  function  tend  to  stabilize  the  instability  in  the  full 
electromagnetic  treatment.  The  theory  suggests  that  prompt  striations  in 
the  magnetosphere  can  form  from  ion  jets  only  if  the  plasma  density  is  suffi¬ 
ciently  tenuous  so  that  electromagnetic  stabilization  is  not  complete. 


C-l.  INTRODUCTION 


Structures  with  magnetic  field  aligned  geometry  have  been  observed 

on  a  time  scale  comparable  to  the  ion-gyrotime  and  have  been  associated  with 

Cl  C2 

plasma  jetting  across  the  ambient  magnetic  field.  *  Called  prompt  stria- 
tions  in  Reference  Cl,  the  linear  analysis  of  these  structures  has  previously 
been  examined  for  an  ion  loss-cone  velocity  distribution  function  in  the 
electrostatic  limit.  5  “  Reference  Cl  showed  that  the  ion  velocity  dis¬ 

tribution  function  attributed  to  the  Buaro  barium  release  was  consistent  with 
the  linear  growth  of  a  purely  growing,  finite  ion-gyroradius  flute  instability. 
Later,  it  was  demonstrated  that  electron-ion  collisions  and  cold  background 

plasma  always  tend  to  stabilize  the  linear  instability,  but  for  the  particular 
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case  of  Buaro,  were  of  insufficient  magnitude  to  provide  complete  stability.  ’ 

In  the  ionosphere  and  magnetosphere,  the  ambient  magnetic  field  strength 
can  vary  over  a  wide  range.  For  example,  in  Buaro  the  ion  kinetic  pressure 
was  much  smaller  than  the  magnetic  field  pressure;  however,  if  a  Buaro-like 
plasma  had  been  injected  into  the  magnetosphere  the  ion  kinetic  pressure  would 
be  larger  than  the  magnetic  field  pressure.  From  a  theoretical  viewpoint  the 
difference  between  the  magnetic  field  strength  in  the  ionosphere  and  magnetos¬ 
phere  implies  that  electromagnetic  effects  on  waves  are  much  more  significant 
in  the  magnetosphere. 

The  purpose  of  this  appendix  is  to  include  full  electromagnetism  into 
the  linear  analysis  of  purely  growing,  finite  ion-gyroradius,  flute  instabili¬ 
ties.  It  is  demonstrated  that  electromagnetic  effects  reduce  the  growth  rate 
from  the  values  calculated  in  the  electrostatic  limit.  Stabilization  by  cold 
background  plasma  and  thermal  spread  in  the  ion  distribution  function  persists 
in  the  electromagnetic  limit.  The  theory  implies  that  a  Buaro-like  plasma 
injected  into  the  magnetosphere  may  lead  to  prompt  striation  formation  but 
only  after  the  plasma  has  become  sufficiently  tenuous. 

Subsequent  portions  of  this  appendix  are  divided  in  the  following  way. 
Section  C-2  gives  the  calculation  model.  Section  C-3  presents  results  from 
the  numerical  evaluation  of  the  electromagnetic  dispersion  relation,  and  Sec¬ 
tion  C-4  describes  application  of  the  theory  to  the  ionosphere  and  magnetos¬ 
phere.  Section  C-5  contains  concluding  remarks. 
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C-2.  MODEL 


As  the  basis  for  the  calculation,  an  infinite  and  homogeneous  quasi¬ 
neutral  plasma  model  is  assumed.  Cold  electrons  and  one  singly  ionized  ion 
species  constitute  the  plasma.  Two  different  ion  velocity  distribution  func¬ 
tions  are  permitted  simultaneously  in  the  plasma.  The  first  assumes  that  ions 
are  cold,  while  the  second  assumes  a  loss  cone  distribution  of  the  following 
form:C1,C4’C5 


-  v  63  /  (v,  -  v,)2\ 

f,(v)  -  g(v.)  h(v  )  * - -V -  exp - - — 2 — )  hCv  )  (C.l) 

1  1  11  2w  Vj  A(6)  \  /  11 

with 

A ( 5 )  =  (|)*  |(1  +  52)  [1  +  erf(6/vS)]  +  (|)*  6  exp(-0.562)J  , 
6=vl/vil  *  (C-2) 


In  Eqs.  (C.l)  and  (C.2),  v^  and  v)(  are  the  perpendicular  and  parallel  velo¬ 
cities,  respectively.  The  parameters  and  v^  are  constants.  Because  the 
present  analysis  deals  with  flute  modes  having  no  variation  along  the  magnetic 
field,  the  dispersion  relation  is  not  sensitive  to  the  precise  form  of  the 
function,  h(v(|),  provided  the  integral  of  h ( v|( )  over  parallel  velocity  is  equal 
to  one.  Hence,  the  function,  h ( v^ ) ,  is  left  unspecified,  and 
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The  wave  electric  fields  are  assumed  to  have  the  plane  wave  form  appro¬ 
priate  for  the  flute  approximation: 


E(x,y,z,t)  =  E  exp[i(kx  -  wt)] 


(C.4) 


In  Eq.  (C.4)  the  x  denotes  the  coordinate  perpendicular  to  the  background  mag¬ 
netic  field  and  parallel  to  the  wavenumber,  k.  The  coordinate  z  lines  in  the 
magnetic  field  direction  and  the  time  is  designated  by  t.  The  wave  frequency 
is  w  and  E  is  a  constant  vector. 

For  the  indicated  plasma  model,  modes  on  the  extraordinary  branch  of 

the  dispersion  relation  with  frequencies  well  below  the  electron-cyclotron 

C6 

of  the  dispersion  relation  satisfy: 
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In  Eqs.  (C.5)  and  (C.6),  the  x-  and  y-directed  electric  field  amplitudes  are 

represented  by  E  and  E  »  respectively.  Subscripts  o,  1,  e,  and  i  denote  cold 

x  y 

background,  fast  ion,  electron  and  ion  quantities,  respectively.  Various  para 
meters  influencing  the  dispersion  relation  are  density  (n),  speed  of  light 
(c),  magnitude  of  electron  charge  (e),  background  magnetic  field  strength  (B), 
and  mass  (m).  Bessel  functions  of  the  first  kind  and  order  p  [i.e.,  J  (b - )] 
and  their  derivatives  with  respect  to  argument  [J  ' (b.)]  play  a  significant 
role  for  the  modes  described  by  the  dispersion  relation. 

The  eigenvalue  condition  is  obtained  by  setting  the  determinant  of 
the  matrix  of  coefficients  in  Eq.  (C.5)  equal  to  zero: 
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For  purely  growing  modes 
speeds  much  smaller  than 
rewritten  as: 


with  linear  growth  rate  y  (i.e.,  u  =  iy)  and  Alfven 
the  speed  of  light,  the  dispersion  relation  can  be 
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In  Eq.  (C.9),  the  symbol  5  is  equal  to  one  for  p=0  and  is  zero  otherwise. 
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In  the  electrostatic  limit  with  cold  background  plasma,  the  dispersion 
relation  depends  upon  only  four  dimensionless  parameters,  p\  x,  5,  and  kb^  = 
kVj/n^.  Equation  (C.8)  indicates  that  electromagnetic  effects  introduce  a 
fifth  dimensionless  quantity,  kc/u^,  into  the  dispersion  relation  and  also 
shows  that  the  electrostatic  approximation  is  valid  for: 
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0  =  EXX  '  kC/“U 


(C.10) 


It  is  useful  to  rewrite  the  dimensionless  quantity  found  in  electromagnetic 
but  not  electrostatic  dispersion  relations  in  the  following  way: 
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Equation  (C.8)  can  now  be  rewritten  as: 


&I  = 


(1+7)?+rt7 


XX 


_  b  c  _ _ 1  1  +  X  _ 

(777  ?  ’  xlxy)'  +  **  + 


+  — 2 — 2  -  U  fd  +  p)  /  +  — ^  -  1  1 

1  +  x‘  XXJ  L  m.  1  +^2  yyj 


(C. 12) 


Equations  (C.10)  and  (C.12)  demonstrate  that  the  electrostatic  approximation 
is  appropriate  in  the  limit  when  the  ion  kinetic  energy  per  unit  volume  is 
much  smaller  than  the  magnetic  pressure  (i.e.,  3^  <<  l). 

Because  it  is  aligned  in  the  wavenumber  direction,  the  E  amplitude 
is  associated  with  an  electrostatic  field  component.  However,  E  is  a  con¬ 
sequence  of  an  inductive  electric  field.  From  Eq.  (C.5)  the  electrostatic 
and  electromagnetic  electric  field  amplitudes  are  related  by  the  expression: 


^x^y  “£x)/exx 


(C.13) 
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In  addition  to  the  x-  and  y-directed  electric  fields,  Faraday's  law  implies 
the  presence  of  a  z-directed  inductive  magnetic  field  with  amplitude 


Bz  =  -i(kc/y)  Ey 


(C. 14) 
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C-3.  RESULTS 


To  determine  the  effect  of  finite  on  purely  growing  finite  ion- 
gyroradius  flute  instabilities,  the  dispersion  relation,  Eq.  (C.12),  has  been 
evaluated  in  terms  of  the  critical  parameters,  x,  kb^,  <$,  and  7.  Results  of 
the  explicit  evaluation  are  summarized  in  Figures  C.l  -  C.4. 

Figures  C.la  and  C.lb  are  plots  of  maximum  values  for  versus  <5  for 

different  values  of  p\  The  numerical  calculations  indicate  that  solutions 

for  positive  3^  are  possible  only  if  e xx  <  0,  which  requires  that  JgCkbj) 

J^kb^)  <  0.  The  difference  between  Figures  C.la  and  C.lb  arises  because  they 

correspond  to  the  first  interval  of  negative  (i.e.,  2.5  <,  kb^  <  3.8)  and 

the  second  interval  of  negative  (i.e.,  5.6  <  kb^  <  7.0),  respectively. 

The  value  of  (B, )mav  depends  on  kb..  For  Figures  C.la  and  C.lb,  it  occurs 
i  max  i 

for  kbj~  3.0  -  3.2  and  6.2  -  6.4,  respectively. 

Both  Figures  C.la  and  C.lb  demonstrate  that  cold  background  plasma 
decreases  the  maximum  beta  permitting  instability.  Specifically,  there  is 
complete  stability  for  the  conditions  given  by  Figure  C.la  for  "p  ,>  .06  and 
in  Figure  C.lb  for  p"  >,  .02.  This  result  is  consistent  with  the  cold  plasma 
stabilization  found  in  the  electrostatic  limit. 

Figures  C.la  and  C.lb  also  show  that  for  specified  values  of  (T  the 
maximum  value  of  3^  which  permits  instability  increases  with  6.  Because 
larger  values  of  6  imply  that  the  loss  cone  ion  distribution  function  is  more 
peaked  and  more  free  energy  is  available  for  instability,  this  result  suggests 
that  finite  3j  decreases  the  tendency  for  instability.  Indeed,  the  numerical 
calculations  show  that  for  fixed  kb^,  6,  and  p\  growth  rate  decreases  as  3^ 
increases.  Maximum  growth  rate  occurs  in  the  electrostatic  limit  with  exx  *  0, 
while  at  marginal  stability  the  maximum  3^  consistent  with  solution  of  the 
dispersion  relation  is  obtained.  This  property  of  the  solutions  is  clearly 
illustrated  in  Figure  C.2  which  is  a  plot  of  maximum  3^  versus  y/n .  for  p"  *  0, 
kbj  *  3.2,  and  for  three  different  values  of  5. 

Thus  far  the  numerical  evaluation  of  the  full  electromagnetic  disper¬ 
sion  relation  has  demonstrated  that  electromagnetic  effects  along  with  cold 
plasma  tend  to  stabilize  purely  growing,  finite  ion-gyroradius,  flute 
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Figure  C.la.  versus  5  for  various  values  of  p 

I  lUa  X 

and  2.5  <,  kjbj  <,  3.8.  (3^^  for  6  • 

is  shown  on  the  right-hand  portion  of 


each  curve, 


w 


instabilities.  However,  the  relationship  between  the  inductive  and  electro¬ 
static  electric  fields  needs  to  be  quantified  in  order  to  ascertain  whether 
or  not  electromagnetic  effects  appreciably  modify  the  wave  electric  field 
polarization.  Mathematically,  Eq.  (C.13)  must  be  evaluated  as  a  function  of  p, 
x,  kbj,  and  <5.  A  typical  result  of  the  calculation  is  presented  in  Figure 

C.3,  which  is  a  plot  of  the  maximum  ratio  of  the  magnitude  of  E  /E  (i.e., 

_  y  x 

lE^Ex!^)  versus  5  for  three  values  of  p  and  2.5  <  kb^  <  3.8.  The  qualita¬ 
tive  features  of  the  instability  for  other  ranges  of  kb^  permitting  exx  <  0 
are  similar  to  the  results  being  presented.  The  figure  shows  that  even  though 
electromagnetic  effects  are  included  in  the  dispersion  relation,  the  wave  elec¬ 
tric  fields  associated  with  the  instability  are  of  a  dominantly  electrostatic 
character  (i.e.,  |Ey/Exlmax  <,0.015).  Moreover,  this  property  of  the  modes 
becomes  more  apparent  as  the  instability  becomes  less  vigorous  either  because 
cold  plasma  density  increases  or  5  decreases. 

Figure  C.4  is  a  plot  of  |Ey'Exlmax  versus  y/n^.  for  6  *  5,  F  *  0.0, 
and  kb^  a  3.2.  The  figure  illustrates  a  typical  feature  of  the  instability, 
i.e.,  peak  values  of  lEy/Exlmax  occur  for  a  growth  rate  intermediate  between 
marginal  instability  and  the  maximum  possible  growth  rate  for  specified  6, 
p,  and  kb^ .  The  physical  reason  for  this  characteristic  of  the  instability 
solutions  is  twofold.  First,  near  marginal  stability  ion  and  electron  "polar¬ 
ization"  currents  are  much  smaller  than  the  electron  EXB  current  (i.e., 

|c  |  >>  | e  | ) .  Second,  at  the  maximum  growth  rate  for  specified  6,  7,  and 

kb,,  the  instability  is  purely  electrostatic  and  e  *  0. 

X  aA 

Having  demonstrated  that  the  instability  electric  fields  are  of  pri¬ 
marily  electrostatic  character,  it  becomes  evident  from  Eq.  (C.14)  that: 

1  «  |BZ/Ey|  =  (aj.j/fi.)  kbjtxf^r1  .  (C.15) 

So  the  most  significant  implication  of  electromagnetic  effects  on  wave  fields 
is  to  permit  the  Bz  field.  Associated  with  this  wave  magnetic  field  is  the 
instability  nature  to  persist  for  a  range  of  finite  3^,  albeit  with  a  reduced 
growth  rate  as  compared  to  the  purely  electrostatic  instability. 


4 

4 
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Figure  C.4.  lEy/Eximax  versus  y/n.  for  kbj  =  3.2, 
6=5,  and  p  =  0.0. 
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C-4.  APPLICATION  OF  THE  THEORY 


Injection  of  plasma  across  the  magnetic  field  in  the  ionosphere  can 

lead  to  the  formation  of  field  aligned  structures  (striations)  with  dimensions 

transverse  to  the  magnetic  field  which  are  smaller  than  the  inverse  ion-gyro- 

radius.  In  particular,  the  Buaro  shaped  barium  release  at  400  kilometers 

altitude  demonstrated  that  these  prompt  striations  can  be  formed  on  a  time 

Cl 

scale  comparable  to  the  barium  ion-gyrotime.  The  electrostatic  approxima- 

5  -3 

tion  is  well  justified  for  the  Buaro  experiment  because  n.  «  4  x  10  cm  , 

_  £  i 

6  =  5,  p  «  .009,  v^  a  1  x  10  cm/sec,  and  B  as  .32  gauss  implies  that  B^»  0.01 
which  is  an  order  of  magnitude  smaller  than  the  maximum  value  of  Sj  ( i . e - ,  0.1), 
permitting  instability  for  6  =  5  and  7  «  0.01  (see  Figure  C -  la ) . 

It  is  instructive  to  calculate  the  minimum  altitude  at  which  complete 
electromagnetic  stabilization  of  the  loss  cone  mode  would  take  place  for  Buaro 
plasma  parameters.  For  this  estimate  it  is  assumed  that  the  earth,  with  an 

3 

approximate  radius  of  6.38  x  10  km,  has  a  pure  dipole  magnetic  field  with 

25  3  r 7 

a  dipole  moment  of  8  x  10  gauss-cm.  Now  Figure  C.la  shows  that  complete 
electromagnetic  stabilization  of  the  flute  mode  for  6=5  occurs  if  Bj  >,  0.2. 

For  the  Buaro  plasma  parameters  a^d  the  model  dipole  field,  this  stabilization 

3 

corresponds  to  B  <,0.08  gauss  c-  equatorial  altitudes  above  3.6  x  10  kilo¬ 
meters.  It  follows  that  shaped  barium  injection  with  plasma  parameters  com¬ 
parable  to  Buaro  will  not  permit  the  formation  of  striations  by  the  finite 
ion-gyroradius,  purely  growing,  flute  instability  in  most  of  the  magnetosphere. 

Although  the  purely  growing  instability  may  not  develop  initially  for 
shaped  barium  injection  in  the  magnetosphere,  it  could  develop  later  as  the 

plasma  density  decreases.  For  example,  consider  shaped  barium  injection  across 

-3  4 

a  magnetic  field  with  strength,  1.6  x  10  gauss,  at  an  altitude  of  3.15  x  10 

C8 

km.  (The  indicated  parameters  are  those  of  an  earlier  magnetospheric  release  .) 

Assuming  that  6=5  and  v.  «  1.0  x  10®  cm/sec,  the  maximum  barium  density  which 

A  2-3 

permits  kinetic  flute  mode  growth  is  n.  a  1.8  x  10  cm  .  For  a  background 

-3  C8  A  p  -3  — 

proton  density  of  16  cm  ,  n^  »  1.8  x  10  cm  implies  p  <  .001  and  stabili¬ 
zation  by  background  plasma  is  ineffective. 
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C-5.  CONCLUDING  REMARKS 


Purely  growing,  finite  ion-gyroradius,  flute  modes  have  been  postulated 

to  be  the  cause  of  prompt  striations  in  a  shaped  barium  release  across  the 

Cl 

magnetic  field  in  the  ionosphere.  Electromagnetic  effects  have  been  included 
in  the  stability  analysis  of  the  flute  modes  in  order  to  ascertain  their  impact 
on  linear  growth.  Although  inductive  electric  fields  are  much  smaller  than 
the  electrostatic  fields  under  all  circumstances,  electromagnetic  effects 
nevertheless  do  modify  the  wave  fields  by  introducing  a  parallel  wave  magnetic 
field  not  present  in  the  electrostatic  limit.  When  electromagnetic  effects 
are  appreciable,  the  ratio  of  ion  inertial  pressure  to  magnetic  field  pressure, 
B^,  is  finite.  Under  such  circumstances,  the  linear  instability  growth  rate 
is  smaller  than  it  would  be  in  the  electrostatic  limit  with  -*•  0.  In  short, 
finite  3^  tends  to  stabilize  the  mode.  As  in  the  electrostatic  limit,  cold 
plasma  and  spread  in  the  ion  loss  cone  distribution  function  also  inhibit  growth. 

Application  of  the  analysis  to  shaped  barium  release  in  the  ionosphere 
and  magnetosphere  has  demonstrated  that  for  a  specific  ion  velocity  distribu¬ 
tion  function  the  total  ion  density  must  be  less  than  a  specific  value  in  order 
that  finite  beta  stabilization  can  be  circumvented.  In  particular,  for  an  ion 

loss  cone  distribution  function  like  the  one  in  the  Buaro  shaped  barium  release, 

5  -3 

densities  less  than  the  Buaro  barium  ion  density  of  4  x  10  cm  are  necessary 
for  the  kinetic  flute  instability  to  be  operative  in  environments  with  B  <,  .08 

3 

gauss  corresponding  to  magnetic  equc:orial  altitudes  ,>  3.6  x  10  km.  Similarly, 

an  example  appropriate  to  the  magnetic  field  strength  and  altitude  of  the  mag- 

re 

netospheric  barium  release  has  indicated  that  barium  ion  r'jnsities  less  than 
2  -3 

1.8  x  10  cm  are  necessary  for  the  kinetic  flute  mode  to  be  unstable.  These 
examples  suggest  that  in  a  weak  magnetic  field  environment,  the  kinetic  flute 
instability  might  be  operative,  not  initially,  but  only  later  as  plasma  becomes 
more  tenuous  and  finite  beta  stabilization  is  reduced.  Moreover,  even  when 
the  instability  is  effective,  finite  ion  inertial  effects  can  reduce  the  growth 
rate  in  a  weak  ambient  magnetic  field  environment. 
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